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Colloidosomes are created in the laboratory from a Pickering emulsion of water droplets in oil. 
The colloidosomes have approximately the same diameter and by choosing (hairy) particles of 
different diameters it is possible to control the particle density on the droplets. The experiment 
is performed at room temperature. The radial distribution function of the assembly of (primary) 
particles on the water droplet is measured in the laboratory and in a computer experiment of a 
fluid model of particles with pairwise interactions on the surface of a sphere. 
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I. INTRODUCTION 

Colloidosomes are hollow spherical structures that are formed by the assembly of colloidal particles at 
the interfaces of two immiscible liquids^. As a result the particles are arranged in a shell that is inherently 
porous. 

The assembly of colloidal particles at liquid interfaces is used in various applications^. Moreover, it is 
a promising technique for the synthesis of novel materials^ and has recently led to the development of 
colloidosomes^, nano-composite particles^, porous solids^, and foams^. 

In this work, we study colloidosomes that are composed of uncharged spherical polystyrene particles 
of /im size moving on the surface of a water droplet in oil. Similar studies have also been done with 
charged particles 7 -. The study of particles on the surface of a sphere dates back to the old Thomson 
problem^ for classical electrons. The statistical physics problem of a one component plasma on a sphere 
has been solved exactly analytically at a special value of the temperature^. Non point-wise particles on 
a sphere have the additional complication of the geometrical frustration, which can be described through 
the so called grain boundary scarsi£~— . There have been attempts to formulate a statistical geometry of 
particle packing^. These systems opened up a new field of research that studies the effect of curvature 
and topology of various surfaces on the organization of matter in a more general sensed. Structuring at 
the surface of a droplet can be viewed as a two dimensional analog to fluid like behavior, crystallization, 
or glass formation in three dimensional systems^. 

Fluids on Riemannian surfaces have been the subject of various studies with few exact analytical 
results^ ' 16 ' 17 , some approximate theories^, and many Monte Carlo (MC) simulations^. 

Colloidosomes with tunable particle density were synthesized experimentally^. A sintering procedure 
is then used to create capsules which can be easily handled. The capsules are then dried to obtain colloidal 
cages. The synthetic details are explained in the next section. 

In this work we give the simplest statistical physics description of the colloidosomc, where we describe 
the interaction of the colloidal particles with the surrounding media, water and oil, simply as a holonomic 
constraint on the particles positions to stay at the water-oil (w/o) interface and treat them as a fluid 
of a fixed number of particles moving on a sphere, the droplet of water in oil, with a mutual pairwise 
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interaction, the pair- potential, at a temperature T. Additional frictional effects have been neglected^. 
The assembly of particles on the sphere is studied both in the laboratory and with a computer experiment 
under certain conditions: number density and temperature. The structural arrangement of the particles 
is characterized through the radial distribution function. The colloidal particles created in our laboratory 
are polystyrene solid spherical hairy particles with controllable diameter of the order of 3 /xm. The 
particles will then exhibit a hard core interaction. Two types of particle pair-potentials were used in the 
Monte Carlo simulation of the fluid, namely the hard-sphere one and the polarizablc hard-sphere one. 

The work is organized as follows: in Section [IT] the colloidosome is described; in Section UTTl the radial 
distribution function as a means to probe the structure of the colloidosome is presented in its mathematical 
definition described in Appendix [XJ its MC estimator, and its experimental measure; in Section [TV] the 
MC simulation results are presented; in Section [V] the theoretically exact results of the MC simulation 
and the experimental results are compared; Section [VII is devoted to concluding remarks. 



II. THE EXPERIMENTAL SYSTEM VERSUS THE STATISTICAL PHYSICS PROBLEM 

The details for the synthesis of the colloidosomes can be found in our previous work^ii. Working at 
room temperature, we first disperse the colloidal particles in a hydrocarbon oil (heptane). Then, water is 
added while the solution is being stirred vigorously. The function of the shear is two fold. It causes the 
water to break up into small water droplets and at the same time it allows to overcome the barrier for 
adsorption of the particles which are assembled randomly at the w/o interface of the droplets. Eventually 
a stable Pickering emulsion^ of water droplets covered by polystyrene (pS) particles, the colloidosomes, 
in oil is formed. The colloidosomes formed have all approximately the same diameter, and this, as well 
as the number of colloidal particles on each colloidosome, does not change after the colloidosomes are 
formed. For further imaging with scanning electron microscopy (SEM) the particles surrounding the 
droplets need to be (partially) sintered in order to form a continuous and stable shell around the droplet, 
a capsule. Heating the mixture to 35 °C for 30 minutes proved sufficient to do so. A small amount of the 
sintered colloidosome dispersion is placed on the SEM sample holder and dried in the fumehood. This 
removes both the water and heptane remnants and leaves only the capsules, the colloidal cages, which 
can then be imaged by SEM. 

The colloidal polystyrene solid particles are synthesized by the dispersion polymerization of styrene 
in alcohol/water—, which is a well established technique for the formation of highly uniform polymer 
particles with a narrow size polydispcrsity. AIBN (azo-bis(isobutyronitrile)) is used as the initiator. 
The presence of a polymeric stabilizing surface functional group is required for a controlled synthesis of 
the particles. A non ionic polymeric stabilizer (poly(N-vinylpyrrolidonc) (pVP)) is used for this study. 
Hence, there is no charge on the surface of the particles. During the polymerization, pVP, which efficiently 
adsorbs on polystyrene, is attached to the particle. In ethanol and water, pVP is soluble. The polymeric 
chains are extended (with a radius of gyration of R g m 15 nm) and are responsible for steric repulsion as 
two particles get in close contact. This guarantees the steric stabilization of the colloidal suspension. In 
the oil phase, pVP is insoluble and the polymer chains are collapsed on the surface of the particle, resulting 
in an attractive potential among the particles. The final dispersion in cthanol/water, therefore, consists 
of polystyrene particles that are sterically stabilized with a layer of pVP. During the polymerization of 
styrene, pVP attaches to the particle by both physical adsorption to the particle's surface and chemical 
grafting. These two mechanisms occur simultaneously, however it is not known to which extent. It is 
believed that the predominant mechanism for stabilization is the physical adsorption. The particles are 
washed with pure ethanol by three centrifugation/redispersion cycles in order to remove the residual 
physically adsorbed pVP. The particles settle due to the centrifugation. The supernatant solution is 
decanted and clean ethanol is added to the remaining particles. The particles are then redispersed and 
the whole procedure is repeated 3 times. The physically adsorbed pVP is removed, which produces a lower 
colloidal stability of the particles in ethanol. Large aggregates were observed during this procedure, which 
is an indication for the presence of an attractive component in the pair-potential between the colloidal 
particles in the suspension. The chemically grafted pVP remains attached to the particle's surface. The 
remaining particles are dried and redispersed in heptane, before the colloidosomes are synthesized. In 
conclusion, the surface chemical properties of the particles arc mainly determined by polystyrene and 
the fraction of pVP that is chemically grafted to the surface, although a precise estimate of the grafting 
density is lacking. It is believed that the grafting density is low, due to the poor colloidal stability of the 
particles in ethanol after removing the physically adsorbed pVP. 

In Pickering emulsions, the particle is adsorbed at the w/o interface and is partly immersed in the 
oil and water phase. The extent of immersion in both phases will eventually have an influence on the 
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particle pair-potential. The three phase contact angle 9 is used to denote the position at the interface 
as shown in Fig. It was determined in our earlier work» and is approximately 9 « 130°, which 

means that the particle is predominantly immersed in the oil. The surface tension of the particle is 
altered by the presence of the surface stabilizing groups, which affects the wetting properties and the 
equilibrium position of the particles at the interface^. An atomistio^ level of description of the core of 
the solid hairy particles suggests the use of Hamaker— calculation for the determination of the interaction 
between the two spherical cores. The calculation predicts an attractive pair-potential which, neglecting 
the detailed behavior close to contact, is proportional to (a/r) 6 , r being the distance between the centers 
of the two cores of diameter a. Unlike this attraction, which is always present, the steric repulsion will 
have a very small range when is obtuse since in this case the particles contact occurs in the oil phased. 
We thus expect the balance between the attractive interaction and the steric effects to depend on the 
angle 9. Moreover other kinds of interactions such as depletion, hydrophobic, solvation, or capillary 
should be taken into account for an accurate description of our system. The simplest description for the 
pair-potential between the particle is the hard-spheres one. 





FIG. 1. The left panel shows a schematic representation of the equilibrium position of the particles at the w/o 
interface; 9 is the three phase contact angle, r is the particle separation, and a is the particle diameter. The right 
panel shows a side view of the particle shell of a colloidosome made with scanning electron microscopy. 



The surface of a sphere of diameter D = 2R is A = ttD 2 . The surface area that a particle, with 
diameter a, can occupy is approximately a = v / 3c 2 /2. The maximum number of particles that can pack 
the surface of a sphere is approximately 



A 2tt 2 



where we assumed that the particles are in a close packing regular hexagonal lattice neglecting curvature 
effects. The maximum reduced particle density on the sphere will then be p m axV < 2/V3 w 1.155. 

Similar experiments^ make use of water and a liquid of higher density, for the initial solution of the 
two immiscible liquids. The droplets in the emulsion will now be of the higher density liquid. In the 
limit of droplets of very high density the particles are expected to be essentially unable to move on the 
droplet. Our working hypothesis will be, instead, to consider the particles as moving freely on the droplet 
surface, completely neglecting the presence of the solvent. We then treat the colloidosome of diameter 
£>, number of particles N, and temperature T, through a canonical ensemble classical statistical physics 
description of the assembly of particles on the water droplet as a fluid of particles constrained to move 
on the surface of a sphere with a pairwise interaction, the pair-potential. 



A. The pair-potential 



Fixing the pair-potential completely defines the fluid model, as described in Appendix [X] and Eq. 
The simplest interaction between two colloidal particles is the hard-spheres (HS) pair-potential 



, m i +oo r < a , . 

^Hs(r) = { Q r > a , (2) 
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where a is the diameter of the spheres and r is the Euclidean center to center distance (see Eq. (J5J). 

The 1 interaction between two neutral particles far apart is dominated by dipolar forces. The simplest 
model potential, suggested by the London forces^, corresponds to hard-spheres of diameter a with 
dispersion attractions, the polarizable hard-spheres (PHS), 

{+00 r < a 

- £ PHs(^) 6 r>a ' (3) 

where epns = ^4^/36 is a positive energy proportional to Hamakcr constant A H which is a property of 
the material of which the particles are made and of the environment where the particles are immersed. 
We here neglect the close to contact — A H a/2A/(r — a) Hamaker behavior as on the average, even in a 
closed packed configuration, there are many more pair interactions between far particles than between 
close ones. 



III. THE RADIAL DISTRIBUTION FUNCTION 



In this work we probed the structure of the colloidosome using the radial distribution function (RDF). 
We compare the experimental RDFs with the ones obtained from MC simulations of a fluid of particles 
moving on a sphere and interacting with a model pair-potential of the kinds described in Section llI A[ This 
procedure will allow us to determine which interaction model best describes the experimental assembly 
of particles. Choosing a as the unit length, the statistical physics problem only depends on the number 
of particles N and the reduced density pa 2 = N/[w(D/a) 2 ] for the athermal HS model and also on the 
reduced temperature fceT/epHS for the PHS one. 



A. In the Monte Carlo simulation 



On a sphere, the Monte Carlo simulation^ solves exactly the statistical physics problem as, since one 
does not have the additional thermodynamic limit problem, it reduces to an integration, as described in 
Appendix iBl 

The particles positions are R = (ri, T2, ■ • • , rjv) with 

Ti = R[sm 9i cos ipii. + sin 9i sin ipiy + cos 0jz] . (4) 
The Euclidean distance between particles i and j is given by 



r 



i j 



Ry/2 - 2r< ■ f j , (5) 



where fj = Vi/R is the versor that from the center of the sphere points towards the center of the ith 
particle. 

The density of particles on the surface of the sphere is 

N 

(6) 



r 4wR 2 ' 

In the MC simulation^ the RDF between two points on the sphere, r and r', is calculated through the 
following "histogram" estimator (see Eq. (|A14[) ) 

5 (d) = ( 5 hist "(d,R)) , (7) 

where d = 2i?sin(arccos(f • f')/2) is the Euclidean distance between r and r', (...) = f„ N exp(— BUnCR,)) 
. . . dR / f c jv exp(— /3[/jv(R)) dR is the thermal average, here 

R 

MR) = 5>(r< 3 -) , (8) 

i<j 

is the total potential energy of the fluid of particles, cj> is the pair-potential, and the integrals are taken in 
such way that £ Sr for i = 1, 2, . . . , N with Sr the sphere of diameter D = 2R, so that dR = Yli dri 
with dri = R 2 dVti = R 2 sin 9id9idipi, and 



„histogram/j r>\ \ " 1 [d-A/2,d+A/2[{nj ) , , 

9 ^ R )-2. N^d) W 
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here lf a ,b[(f) = 1 if r e [a, b[ and otherwise, and 



( 



d + A/2 
2~R 



) 



2 



( 



d- A/2 
2i? 



) 



2 



n ld (d) = TV 



(10) 



is the average number of particles on the surface [d — A/2, d + A/2[ for the ideal gas of density p. p 2 g(d) 
gives the probability that sitting on a particle at r one has to find another particle at r'. 



B. In the experiment 

The positional data of the particles in the colloidal cages is directly extracted from SEM images, which 
allowed the calculation of the particle separation for all visible particle pairs. A is set to an arbitrary 
value of a/20. To exclude edge effects, a selection of particles located sufficiently at the center of the 
SEM image of the colloidal cage is taken into account. The RDF is determined from just one hemisphere. 
The particle positions from 5 SEM images of similar colloidal cages were used for the statistical average. 
The detailed procedure, the selection of particles, and validation of the procedure is described in our 
previous work>2£. In that work we calculated the radial distribution function from just one SEM image. 
Here we refined that analysis averaging the results from 5 SEM images which is in spirit closer to the 
procedure used in the MC simulations. Although five images are still a rather small number our present 
procedure carry nevertheless more information than the one used in paper—. The absolute error on the 
experimental g(r) is around 0.3. In the experimnt unlike in the simulations each image measurememnt 
is scorrelated from the other. 



IV. MONTE CARLO RESULTS 

Wc performed constant N , p, and T canonical MC simulations^. A typical run would consist of about 
5 x 10 5 iV single particle moves, keeping the acceptance ratios constant (w 0.3). In all the presented 
graphs of the simulated RDF the statistical error from the MC integration are not visible on the chosen 
scale. 

We initially chose the PHS model pair-potential to see how the RDF would change upon changing 
the temperature and the density. In order to find agreement with the experimental results it proved 
necessary to use the simpler HS model, as shown in Section [V] We then compared the HS results with the 
soft-sphere model </>ss( r ) = ess (7) considered in Ref. [l|| For case "a" of Table U a reduced temperature 
of fcsT/ess = 0.05 is sufficient to have similar RDFs for the HS and the SS model on the half hemisphere, 
but when looking at the RDFs on the whole sphere, the SS RDF, unlike the HS one, shows relevant 
correlations between particles at opposite poles (\g(2R) — 1| w 0.3). 



A. Dependence on the temperature 

The HS model is athermal so the structure is independent of temperature but only depends on the 
density. We thus simulated the colloidosome "c" in Table U with the PHS model. And we chose different 
values of the reduced temperature, fcsT/epHS, to see how the RDF would change. 

As expected we found the occurrence of an ordered structure at small reduced temperatures (see Fig. 
[2]) . In particular we observe the formation of a regular hexagonal lattice distorted by the curvature of the 
spherical surface. Recall that in a planar perfect hexagonal lattice arrangement of spheres of diameter a 
the first coordination shells are as follows: r/cr = 1, 1/3, 2, y/7, 3, 2^, Vl3, 4, Vl9, a/21, 5, . . .. 
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FIG. 2. (Color online) RDF for approximately one hemisphere of the colloidosome "c" of Table|T]with N = 1498, 
D/a — 23.8 and various reduced temperatures, calculated with MC simulations of the PHS fluid. Also shown are 
the locations (vertical lines) of the first eight coordination shells of a regular planar hexagonal lattice of the hard 
core particles (here a — 2.72 (jm). The mismatch between the peaks of the RDF and these shells is a manifestation 
of the curvature of the surface. 



From Fig. [2]we can clearly see how at this reduced density, 0.84, well below the maximum density, the 
PHS model reduces to the HS model for reduced temperatures > 1. As we lower the temperature, the 
attractive tail in the pair-potential starts to play a role resulting in a solidification of the fluid. As the 
fluid crystallizes, it may experience the cage effect going through glassy phases. The particles become 
confined in transient cages formed by their neighbors. This prevents them from diffusing freely on the 
surface of the sphered. A related problem is the extremely long MC equilibration time necessary to draw 
the RDF of the figure at a temperature of 0.1, starting from a disordered initial configuration. In the 
limit of T — >• 0, in our calculation, the equilibrium configuration of the (classical) particles is the one Ro 
for which [7/v(R) has its minimum: the probability density is zero everywhere except on Ro. We have a 
spontaneous breaking of the rotational symmetry (see Appendix [S| . The monotonously increasing tails 
in the PHS pair-potential produce an equilibrium configuration with the particles forming one cluster of 
touching spheres. On the contrary in the SS model the equilibrium configuration will be one where the 
inter-particle spacing depends on the density. 



B. Dependence on the density 



For case "c" in Table Q] (D/a = 23.8 and a = 2.72 /im) we chose different values of the density to see 
how the RDF would change for the PHS model at a relatively high value of the reduced temperature 
^s^/epHs = 9.1. 

We succeeded in reaching high particle densities (without overlaps) by placing one particle at the 
north pole and then others centered at 8 = 2narcsin(l/2i?) and (p = 2marcsin(l/2i?sin#) with n,m = 
1, 2,3,.. .. This way we were, in particular, able to reach the 0.91 critical density observed by Prcstipino 
Giarritta et alr^- for HS. In doing so we observed the splitting of the second peak into a pair of adjacent 
peaks corresponding to the second and third coordination shells of a regular hexagonal lattice (see Fig. 
ED- 
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FIG. 3. (Color online) RDF for approximately one hemisphere of the colloidosome with D/a = 23.8 and a — 
2.72 /im at a reduced temperature fcsT/epHS = 9.1 and various densities, calculated with a MC simulation of the 



From Fig. [3] we can clearly see how the fluid tends to reach an ordered phase at high densities (even 
at high temperatures). The realization of these ordered phases will go through the formation of colloidal 
geometrical cages (due to geometrical frustration) on the surface of the water droplet which is inevitable 
as the density slowly approaches the maximum density at any temperature. 



V. COMPARISON BETWEEN EXPERIMENTAL AND MONTE CARLO SIMULATION RESULTS 

The results of the experimental colloidosomes are now compared with the Monte Carlo (MC) simu- 
lations. Scanning Electron Microscopy (SEM) images of two different colloidal cages can be found in 
Fig. 2J The RDF of these colloidal cages and two others are shown in Figs. [5]and|ni The experimental 
colloidosomes studied differ one from the other by particle size and particle density, the water droplets 
were of the same diameter and the temperature was room temperature, as summarized in Table HI The 
same values for number of particles, N, and sphere diameter, D/a, are used in the MC simulations. Our 
first choice for the pair-potential was the HS fluid model, as justified in Section [TT1 

It is important to stress that in the experiment we measured the RDF from 5 images of different colloidal 
cages. Now, there are two processes responsible for the assembly of the particles on the colloidosome: i. 
the adsorption of the particles on the interface at the moment of the formation of the Pickering emulsion 
and ii. the motion of the particles on the interface. Our experimental measure is clearly not able to 
discriminate which one of the two processes is the more relevant, even if we expect that the structure 
of the colloidal cages obtained after the sintering procedure will carry no history of the former process. 
Moreover theoretical studies of process i. are, to the best of our knowledge^, much less developed than 
the ones of the latter. Our computer experiment only takes into account the second process assuming 
that the colloidosome is formed and the particles are in thermal equilibrium on the droplet. 

From Figs. [5] and [6] we can see that there is a good agreement between the experimental and the 
theoretical RDF. This indicates that the HS fluid model gives a good description of the experimental 
system. The snapshots from the MC simulation of the colloidosome differ from the SEM images of the 
colloidal cages. The colloidal cages are formed by a network of touching particles. The structure of 
the experimental fluid points to a pretty strong short-range attraction between the particles mainly as a 
result of the sintering process. The measure of 'structure' used, g(r), is not sensitive to these structural 
differences. 



PHS fluid. 





case a (/im) N D/a pa 2 

a 4.80 561 13.5 0.98 

b 3.32 1065 19.5 0.89 

c 2.72 1498 23.8 0.84 

d 2.56 1449 25.3 0.72 



TABLE I. Characteristics of the experimental colloidosomes analyzed. In all cases, the water droplet was of the 
same diameter D = 64.8 fim. Different colloidosomes differed by the diameter a of the colloidal particles and 
by the number N of colloidal particles they carried. The same systems, "a" , "b" , "c" , "d" , have been studied 
through MC simulations. 



We also simulated the experimental colloidosomes with the more realistic PHS model. Initially, the fluid 
with the highest particle density ("a") was used to adjust the reduced temperature. By trial and error we 
found that fc^T/epHs = 0.3 gave a satisfying agreement with the experiment (see Fig. [7]). However when 
we simulated the other colloidal cages with the same reduced temperature (the experimental temperature 
in all cases did not vary and the Hamaker constant did not change from one colloidosomc to the other) 
we found disagreement between the MC simulation and the experiment as clearly shown by the last panel 
of Fig. [7] This is an indication that the particles used in the experiment do not interact as PHS. An 
explanation for this is the balance, in the oil phase, between the steric repulsion of the polymer chains 
and the Hamaker attraction. However, steric repulsion through the oil phase is unlikely, because pVP is 
insoluble in heptane. Another, possible explanation is that during emulsification the attractive interaction 
is balanced by the shear that is applied, and this could be reflected on the capsules structures sintered 
after the emulsification process. 

In Table [H] we report the excess internal energy per particle measured in the MC simulations of the 
PHS model in the various systems studied. We can clearly see that as the fluid develops towards a solid 
phase there is a lowering of the energy. 
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FIG. 5. (Color online) The experimental (dashed red line) and simulated (solid black line) RDF of colloidosomes 
"a" (top panel) and "b" (bottom panel) of Table [I] The fluid model used in the MC simulations was the HS. 
The graphical representation of a snapshot of the particle positions when the MC has reached equilibrium, shows 
resemblance with the SEM images. 
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FIG. 6. (Color online) Same as Fig. [5] for the colloidosomes "c" (top panel) and "d" (bottom panel) of Table [I] 
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FIG. 7. (Color online) The experimental (dashed red line) and simulated (solid black line) RDF of colloidosomes 
"a" , "b" (top panel) and "c" , "d" (bottom panel) of Table [I] The fluid model used in the MC simulations was 
the PHS with a reduced temperature of 0.3. 



fcsT/epHS 


N 


D/a 


u ex /cphs 


0.3 


561 


13.5 


-2.1509(6) 


0.3 


1065 


19.5 


-1.862(3) 


0.3 


1498 


23.8 


-1.732(2) 


0.3 


1449 


25.3 


-1.586(3) 


0.1 


1498 


23.8 


-2.3484(8) 


1.0 


1498 


23.8 


-1.5450(4) 


9.1 


1498 


23.8 


-1.5136(3) 


9.1 


1747 


23.8 


-2.0932(2) 


9.1 


1248 


23.8 


-1.0741(2) 



TABLE II. Excess internal energy per particle u ex — {Un)/N for the simulated PHS fluids on the sphere. 



VI. CONCLUSIONS 

We have studied a colloidosome of polystyrene hairy particles of size moving on the surface of a 
water droplet in oil, both experimentally and theoretically through canonical Monte Carlo simulations, 
which is the ensemble of choice for the description of the experimental system, where the number of 
particles does not vary. In particular wc studied the radial distribution function. Agreement was found 
between the experimental measure and the measure in the computer experiment of the theoretical model 
of the fluid of pure hard-spheres particles on the surface of a sphere. We did not find agreement between 
the experiment and the polarizable hard-sphere fluid model. 
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It would certainly be interesting to pursue a different determination of the radial distribution function 
through the imaging of the same colloidosomc without going through the sintering procedure. This would 
allow an unbiased determination of the structure of the fluid of particles in their thermal equilibrium on 
the droplet surface. 

Within the Monte Carlo simulation, a wide range of particle densities on the colloidosome was studied. 
At high density, the particles tend to arrange in a hexagonal lattice, distorted by the curvature of the 
droplet and the radial distribution function shows clear signatures of the first coordination shells. While 
at low densities a fluid like behavior is manifested. 

Our Monte Carlo simulation results further show that the addition of an attractive tail to the pure 
hard-sphere pair-potential allows to reach solidification by lowering the temperature even at low densities. 
We discussed that for pair-potentials with support on the whole [0, 2R], crystallization has to be expected 
at low temperatures at any density. 

In our Monte Carlo study we have only considered direct interactions between the colloidal particles 
and not solvent mediated indirect interactions like excluded volume depletion forces, the Duits interaction 
between hairy hard-spheres^, or the Gourney solvation interaction, which depend on the thermodynamic 
state of the system^. We have simulated two fluid models on the spherical surface: the athermal HS one 
and the PHS one. By tuning the reduced temperature in the PHS model so as to get a structure similar 
to the one of the experimental case "a" of Table HI we were then unable to reproduce the experimental 
radial distribution function of the other cases "b" , "c" , and "d" . Only the HS model agrees with all four 
experimental cases. 

Our simulations shows that the HS fluid model has small correlations between particles at opposite 
poles even at high densities^, this is not anymore so for particles with a soft corei^. From the point 
of view of our work this remains just a theoretical prediction as our experimental measure of the radial 
distribution function is only able to probe half hemisphere. Moreover we expect the polystyrene particles 
used to be well described by the hard core pair interaction. A further interesting comparison between 
the experiment and the simulation would be to compute the orientational correlation function Q(r^. 

A possible further development of the work could be the realization of the binary mixture of small and 
large particles on the water droplet^ to find experimental evidences for demixing predicted by the non 
additive hard sphere model with negative non additivity^. Or as a possible way to push the fluid-like 
behavior at larger densities, diminishing the glass gap, as predicted by the additive hard sphere model^i. 
A natural extension would then be the multicomponent mixture which in its polydisperse limit may leave 
no space to the glass. It would also be possible to simulate the particles as penetrable-square- well ones^. 

Colloidosomes may be used to isolate viruses when at the moment of the formation of the Pickering 
emulsion only one is captured inside each droplet. This may be a way to overcome the usual staining 
procedure. Of course the opposite situation may also be possible when many living cells, like for example 
eukaryotic flagellar, coordinate themselves in the confined geometry of the drop. 



Simplicio: As for me, I understand in a general way how the two kinds of natural motions give rise to the 
circles and spheres; and yet as to the production of circles by accelerated motion and its proof, I am not entirely 
clear; but the fact that one can take the origin of motion either at the innermost center or at the very top of 
the sphere leads one to think that there may be some great mystery hidden in these true and wonderful results, 
a mystery related to the creation of the universe (which is said to be spherical in shape), and related also to the 
seat of the first cause. 

Salviati: I have no hesitation in agreeing with you. But profound considerations of this kind belong to a 
higher science than ours. We must be satisfied to belong to that class of less worthy workmen who procure from 
the quarry the marble out of which, later, the gifted sculptor produces those masterpieces which are hidden in this 
rough and shapeless exterior. Now, if you please, let us proceed. 

[■■■] 

Sagredo: I really believe that just as, for instance, the few properties of the circle proven by Euclid in the 
Third Book of his Elements lead to many others more recondite, so the principles which are set forth in this little 
treatise will, when taken up by speculative minds, lead to many another more remarkable result; and it is to be 
believed that it will be so on account of the nobility of the subject, which is superior to any other in nature. 
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Appendix A: The pair correlation function 



Given a classical system of N particles of mass m moving in a region 1Z of a Riemannian manifold of 
dimension d and metric tensor g MI /(q) with Hamiltonian, 

Hn = Tn+U Ni (Al) 
1 N 

Tn= 2^Y,9^ (<K )JVPi" . ( A2 ) 

i=l 

U N = ^C/jv(qi,...,qjv) , (A3) 

where we denote with q = (g 1 , . . . , g d ) a point of the manifold, with q^ = (qj, . . . , qf) the coordinates of 
particle i and with pi = (pu, . . . ,Pid) its canonically conjugate momenta, and we use the Einstein sum- 
mation convention to omit the sum over the repeated Greek indices. The canonical ensemble probability 
density to find the statistical system of distinguishable particles, the fluid, in thermal equilibrium at an 
inverse temperature /3 = 1/ksT (with ks Boltzmann constant) with coordinates Q = (qi, . . . , qjv) and 
momenta P = (pi , . . . , pjv) is 

^ (Q ' P ^ ) = ^^4n e ^ (rjV+ " N) ' (A4) 

where h is Planck constant and the normalization factor is the partition function of the canonical 
ensemble of the identical particles 

where in the second equality we performed the Gaussian integral over the conjugated momenta so that 
dV = dvi • • • dvN with dvi = ^11^=1 (here g = det ||<7^(q;)|| = [det 1 1^^" (qi) 1 1] —1 ) the infinitesimal 
"volume" element of the manifold and A = y/2ir/3h 2 /m is the de Broglie thermal wavelength. To justify a 
classical treatment of the statistical properties it is necessary that A be much less than the mean nearest 
neighbor distance between the particles. In the last equality we used the definition of entropy and F is 
the Hclmholtz free energy. 

If the particles move on a sphere of radius one, 1Z = S±, then the coordinates are the polar coordinates 
on the sphere q^ = tfi) with 9i € [0, tt] the polar angle and ipi € [0, 2tt) the azimuthal angle, and the 
metric tensor is 

||<fc*(q)l| = ( I °2 J > (A6) 
V sm J 

so that det ||(7 Mt /(qi)|| = sin 2 9i and dv; = sin9id9id(fi = ddi the solid angle spanned by the three 
dimensional vector of the position of particle i in the Euclidean space with origin on the sphere center. 
Given a physical observable only function of the coordinates On = Ojv(Q), we can then measure its 
average value as 



For example the one body correlation function^ for the particles on a sphere of radius R is measured 



as 

/ N 



where J( 2 )(q, q') = 6(9 — 9')8{tp — ip')/y/g is the Dirac delta function on the manifold. We now use the 
fact that our potential energy ([5} is invariant under any rotation of the reference frame to say that p(q) 
has to be independent of q and (by integrating (|A8I) over civ) we must have p(q) = p = N/ (AttR 2 ). 
The two body correlation function^ is measured as 

ff(q,q) = (2^ — ^2 R 2 ) I [p(q)p(q )] ■ ( A9 ) 

\ i=£j I 
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Because of rotational invariance g can only depend on the geodesic distance d between q and q'. We can 
then calculate it on a reference frame where ip = <p' so that d = R(6 — 6') and 

9( d ) = (E %~ Bi } S{ff r^ S^ - - W )\ /(i?V) • (AlO) 

\ /t^ sin 6/j sin / ' 

If we now choose 6 = 6 + 6' and integrate over dv' we get 

where 0^ = 6i — 8j and ^ = tpi — ipj . We can use rotational invariance to choose the sphere north pole 
sitting on particle j to get further 

g{d) = (^ 8(6 - OMipij)^ /(NpR 2 sm6) . (A12) 

In place of the geodesic distance R9 we can use the Euclidean distance d = 2Rsm(6/2). We can then use 
the equality 5(6 — 6^) = 5(d— dy)i£eos(0/2), here dij = 2i?sin(%/2), to write 

g(d) = (^(d-diMVii)) /{Npd) 

= r£5(d- riJ )5{^)\ /(Npd) , (A13) 

where r$j is defined in Eq. ([5]) of the main text. Now we use rotational invariance noticing again that 
given any two point q^ and q^ on the sphere one can always find a reference frame in which ifi = (fj to 
get 

g(d) = / £ 8(d - r«)\ / (Np2nd) . (A14) 



Appendix B: The Monte Carlo simulation 

In the Monte Carlo integration one does a random walk^i in Q with 6i G [0, tt], <pi E [—tt/2, 7t/2) for all 
i = 1, . . . ,N with periodic boundary conditions: <p = ip + 2ir and 6 = 6 + tt. In the Metropolis algorithm^i 
one takes as the acceptance probability 

A[Q ^ Q'] = min { 1, e -^[^(Q')-^(Q)] H|i sin ^ 1 (B1) 
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